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Abstract 

An interesting description of a collinear configuration of four par- 
ticles is found in terms of two spherical coordinates. An algorithm 
to compute the four coordinates of particles of a collinear Four-Body 
central configuration is presented by using an orthocentric tetrahe- 
dron, which edge lengths are function of given masses. Each mass 
is placed at the corresponding vertex of the tetrahedron. The center 
of mass (and orthocenter) of the tetrahedron is at the origin of co- 
ordinates. The initial position of the tetrahedron is placed with two 
pairs of vertices each in a coordinate plan, the lines joining any pair of 
them parallel to a coordinate axis, the center of masses of each and the 
center of mass of the four on one coordinate axis. From this original 
position the tetrahedron is rotated by two angles around the center 
of mass until the direction of configuration coincides with one axis of 
coordinates. The four coordinates of the vertices of the tetrahedron 
along this direction determine the central configuration by finding the 
two angles corresponding to it. The twelve possible configurations 
predicted by Moulton's theorem are computed for a particular mass 
choice. 
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1 Introduction 



The main purpose of this paper is the efficient computation of the 
colhnear central configurations corresponding to a given choice of four 
masses. A general introduction to the subject of central configurations 
is the Saari's book [I]. The collinear four body central configurations 
is a particular case of the general collinear N-Body central configu- 
rations studied by Moulton [2]. He determines one half of A^! classes 
of different collinear configurations of N particles, and a formal con- 
struction of each by a long-time consuming iterated algorithm. The 
only three-dimensional Four-Body central configuration is an equilat- 
eral tetrahedron found by Lehmann-Filhes [3]. An algorithm to com- 
pute the non-collinear Four-Body planar central configurations was 
distributed in the web |4J. The collinear Four-Body central configu- 
rations are now considered in this paper using again the coordinate 
system published in reference [5]. This set of coordinates seem to be 
very useful for the study of these problems. 

We have computed previously central configurations of the Four- 
Body Problem for the three-dimensional and the planar non-collinear 
cases in [5] , [6] , and [1] . It is natural to consider at present the collinear 
case that was previously ignored. The new four-body coordinates is 
again the tool to attain the task [5]. That coordinate system for 
collinear problems is so much simplified that the first action is to 
present that reduction in the number of coordinates. 

The masses of the four particles mi, m2, m-a and m4 are positive, 
generally different, but the values could be repeated. An Euclidean 
inertial system of coordinates is the frame to determine the position 
coordinates of the four particles with the center of mass of them at 
the origin of coordinates. Then with no loss of generality the four 
particle positions are function of nine generalized coordinates, which 
for the collinear case reduce to only four coordinates. Starting from the 
general case, the Euclidean coordinate system is transformed to the 
frame of principal axes of inertia [7j by means of a three dimensional 
rotation G parameterized by three independent coordinates. Three 
more coordinates are introduced, as scale factors Ri, R2, R3, which are 
three directed distances closely related to the three principal inertia 
moments through 

h=fiiRl + Rl), l2 = ti{Rl + Rl), and /g = + ^i) , (1) 
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where /x is the mass 
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With the first rotation and the change of scale the resulting four- 
body configuration has a moment of inertia tensor with the three 
principal moments of inertia equal. A second rotation G' does not 
change this property. 

The cartesian coordinates of the four particles, with the center of 
gravity at the origin, written in terms of the new coordinates are 
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where G and G' are two rotation matrices, each one a function of 
three independent coordinates such as the Euler angles, and where 
the column elements of the constant matrix 



E 



^ ai a2 as 04 
h 62 ^3 64 

\ Ci C2 C3 C4 



(4) 



are the coordinates of the four vertices of a rigid orthocentric tetrahe- 
dron, with the center of mass at the origin of coordinates, namely: 



aiuii + 02771-2 + 03777,3 + 047774 = , 
6i777i + 62"i2 + bsms + 647774 = , 

ClTTll -I- C27772 + 637773 -|- C47774 = . 

We used the following notation for the matrix 
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An equivalent condition in order to have three equal inertia mo- 
ments for the rigid tetrahedron is expressed as 
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An orthocentric tetrahedron has the property that the perpendic- 
ular Hnes to the faces trough the four vertices intersect at the same 
point. Orthocentric tetrahedra were considered by Lagrange in 1773 
[S] . Other old references on orthocentric tetrahedra are found in a pa- 
per by Court [9], he calls them orthocentric and orthogonal. Placing 
the four masses at the corresponding vertices, the intersection point 
is actually the center of mass of the four masses, and the moment 
of inertia tensor of the four particles has the same principal value in 
any direction. Equation ([7]) imply that the inertia tensor of the rigid 
tetrahedron is proportional by a factor 2// to the unit matrix. 

The proof of these properties is found in reference [5], where we use 
the notation m = mi + m2 + + 1114 for the total mass of the system. 
It was shown that the position vector of one vertex is orthogonal to 
the three vectors between two vertices of the opposite face (to the 
position vector of the vertex.) The square of the distance between 
two vertices is given by 



This last is the condition to have a moment of inertia tensor with 
the same three principal moments of inertia. The six edges of the 
tetrahedron should be equal (for an arbitrary to the square root of 
the right-hand side of this equation. The volume of this tetrahedron 
is equal to 1/6 if ^ is selected as above. 

A different important coordinate system to fix the origin for mea- 
suring the G' rotation is presented in this preprint. Two other coor- 
dinate systems were presented in references [4| and [5j. In this paper 
the original position of the orthocentric tetrahedron is chosen with 
two pairs of particles each in a coordinate plan, the line joining two 
of them orthogonal to a coordinate axis, the center of masses of each 
and the center of mass of the four coinciding with one coordinate axis. 
More explicitly particle 1 has the coordinates 






(9) 



Particle 2 has the coordinates 




(10) 
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Particle 3 has the coordinates 



(«3,&3,C3) 



Particle 4 has the coordinates 
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(11) 
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This rigid tetrahedron is the generalization of the rigid triangle of 
the Three-Body problem with the center of mass at the orthocenter 
discussed previously in [10] . The same triangle was used with different 
purposes by C. Simo [11]. 

Assuming the four particles move on a straight line, Moulton [2] 
finds the motion of the line is in a constant plane. The number of 
Euclidean coordinates changes from twelve to eight. The first rotation 
around the center of mass is a function of one angle ip that rotate the 
line from its position in the Euclidean plane to the one dimensional 
line where the relative positions move. The three scales Ri are reduced 
to one scale Ri = R. The other two scales of the general system are 
zero: R2 = R3 = 0. The matrix of scales become proportional to 
the projector along the direction (1,0,0). The second rotation G' lost 
the angle around the last straight line direction, and therefore it is 
a function of two angles 9,<j), the angles which are obtained rotating 
from the original position of the tetrahedron to the position where 
the third coordinate coincides with the scaled coordinate along the 
straight line, namely 

/ Xi X2 X3 X4 

yi y2 2/3 2/4 
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where G' is parameterized by 
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This is also expressed as 



xi X2 X3 X4 \ I Rcostp \ 

Vi y2 Us \ = \ R sin ip ( cos 6 sin 9 cos (p sin 6 sin "j E . 

yooooy \ J 

(15) 

The number of independent coordinates for a coUinear configura- 
tion is four. The coordinate V determines the rotation of the hne. 
The coordinate R measure the expansion of the particles in the line. 
Its inertia moment is equal to fiR^. The geometric configuration is 
determined by the two angles 6 and (p 

^ ri r2 ra r4 ^ = ^ cos 6* sin0cos(/> smOsmcf) ^E, (16) 

equivalently 

Tj = Oj cos + hi sin 9 cos </> + q sin 9 sin . (17) 

These are the coordinates along the collinear direction where the pro- 
jection of the orthocentric tetrahedron is performed. They obey the 
normalization 

4 

ruirf = fi , (18) 

1=1 

which is a consequence of equation ([7]) . 

The hemisphere of the two coordinates 9 and (p is sectioned in 
twelve disjoint spherical triangles separated by the six main circles in 
the intersection of the sphere with the planes touching the orthocen- 
ter and orthogonal to the edges of the orthocentric tetrahedron. Each 
of these circles is the geometric set where two coordinates rj have the 
same value. On both sides of these circles the two corresponding coor- 
dinates obey a different inequality. Each spherical triangle is the set of 
points with a different permutation of the size of the four coordinates 
Tj. See figure 1 to illustrate these properties. 



2 Equations of Motion 

I assume for simplicity that the potential energy is given by the New- 
ton potential (the gravitational constant is equal to 1) 

K=-ii:^, (19) 

R ~r \ri — rj\ 
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r(7) > r(13 

r(6)>r(13) 

Figure 1: The set of the four-body coUinear configurations in the central circle 
by the stereographic projection of the hemisphere of the two angles position 
of the orthocentric tetrahedron. The great circles represent the r.j positions 
where two particles have the same coordinate. Each point on the hemisphere 
represent a different coUinear configuration. The particular permutation is 
determined by the inequalities of the coordinates associated to two masses 
on both sides of the great circles. Each of the interior of the twelve spherical 
triangles represent the set of different coUinear configurations with the same 
permutation order. The values of the masses are mi = 20, m2 = 13, 7713 = 
7, m4 = 6. The inequalities between coordinates in figure are labeled with 
the numerical value of the corresponding masses. To simplify compilation 
the r's in figure are not italic fonts. 
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although our results may be generalized for any potential with a given 
power law of the relative distances between particles R |rj — rj\. 

I also computed the kinetic energy as a function of the new coor- 
dinates, which is given by 



R' + R^dj"^ + ni + ni] 
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R^ + R^i^p^ + <j)'^ sin^e + e^) 

(20) 

where ip is the angular velocity of the first rotation, and (Oi, ri2, ^^3) = 
((/> cos 9, —(f) sin 9,0) is the corresponding angular velocity vector of the 
second rotation. 

The equations of motion follow from the Lagrange equations de- 
rived from the Lagrangian T — V as presented in any standard text 
on Mechanics [7], |12) . 

However, the tp coordinate related to the first rotation produces 
conservation of the angular momentum vector in the inertial system 



dT 

—J = flR^Tp = 

oyj 



(21) 



where is the numerical value of the conserved angular momentum. 
The Lagrangian equation of motion for the scale coordinate is 



d'^ ■ ■ ■ rIV 

iJi^R + ^,R{i,^ + 02 sin2 + 9^) = -^, 



(22) 



The Lagrangian equations of motion for the two coordinates asso- 
ciated with the second rotation are 



(23) 



and 



^ifiR^) = fiR^'^ sm9cos9-^ 
dt^ ^ d9 



(24) 



There is one more constant of motion, namely the total energy 



T + V = l 



R^ + R^{i)'^ + (t)^ sin^0 + 



+ V. 



(25) 



Writing the kinetic energy in terms of the angular momentum con- 
stant of motion instead of the ip velocity lead us to 
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The energy conservation is thus expressed as 



2 
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where V represents the potential energy. 



3 Central configurations 

The colhnear central configurations are characterized in our coordi- 
nates by constant values of the 9 and (j) angles. For these cases the 
angular velocity vector Vt is the null vector, and the conservation equa- 
tions of moment and energy, (pT|) and (p7|) respectively, become 



Pw, = ^-R^iJ . (28) 



and 



Note that the potential function V becomes equal to a constant di- 
vided by coordinate R. 

These equations are identical to similar equations obtained for the 
Euler and Lagrange central configurations of the Three-Body problem 
|13j . They are formally the same as the equations for the conies in the 
Two-Body problem in terms of the radius R and the true anomaly ip. 

The constant values of the angles 9 and ^ referred to above are 
not arbitrary but they are determined by the two equations of motion 
([231) and ^ 

dV dV 

These two equations are equivalent to the equations for central 
configurations that are presented generally as Moulton did [2j , written 
in terms of our notation as the four equations 

dV 

— = Xrnin {i = 1, 2, 3, 4) , (31) 

where A is the same parameter for the four equations. Using that the 
sum over the index i is zero on both sides of this equation because the 
origin of coordinates is the center of mass, it follows that only three 
of them are independent. Canceling A one has only two independent 
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CoUinear central configurations with masses 20, 13, 7, 6. 
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= 0.246199851994839 


r(13) = - 


-0.288938654946377 


r(20) = 


-0.189462920315934 


r(6) = - 
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r(6) = 
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= 4.06952228046042 


r(20) = 


0.372886439009567 


r(7) = 


= 0.602862162854886 


r(7) = 


0.022933465266316 


r(13) = 


= 0.259974993861172 


r(6) = - 


-0.200419706448564 


r(6) = 


-0.027293561159533 


r(13) = - 


-0.493518830643398 


r(20) = 


-0.371797434661112 


<P 


0.990022337434498 
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4> 


= 1.38056997484206 
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r(7) = 
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= 0.619712718177074 


r(20) = 


0.166846782575173 


r(6) 


= 0.35795101360832 


r(6) = - 


-0.178777559500521 


r(13) = 


= 0.064367025252091 


r(13) 


= -0.486911083794 
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0.011823000179912 


r(20) = 


= 0.162048815317523 
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Figure 2: Stereographic projection of the hemisphere of the two angles motion 
of the orthocentric tetrahedron. The great circles are the positions where two 
particles have the same coordinate. Double intersection with a right angle 
of these circles are double-double coUisions. Triple intersection is at triple 
collision. Each of the twelve spherical triangles have a right angle and two 
vertices of triple collision. The values of the masses were mi = 20, m2 = 
13, ms = 7, m4 = 6. The centers of the small circles are at the positions 
where the central configurations were computed in the table. 
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Figure 3: Two different times for the motion following a four-body si- 
multaneous collinear central configurations. The values of the masses are 
mi = 20, m2 = 13, ma = 7, m4 = 6 with the permutation order 13, 7, 20, 6. 
The ellipses are labeled with the numerical value of the corresponding masses, 
eccentricity was chosen 0.7 to enhance the elliptic shape. The particular per- 
mutation from the table of central configurations is chosen disregarding the 
cases where one ellipse is very small, and the cases when two ellipses are very 
same equal, as is deduced from the values of the tabulated coordinates r,. 
The center of mass is at the common focus, at the intersection of the two 
straight lines. 
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equations. The two independent equations are obtained from equation 
(jlSp by taking the derivative of both sides with respect to 9 and (j) 



^ dr ^ 
1=1 1=1 




(32) 



Substitution of equation ([3T]) in these equations lead to the equivalent 
conditions ([30]) . 

In order to compute four body collinear configurations one uses 
those conditions (I30p by defining the two 4- vectors Si = (dri/dO) and 
sinOti = [dri/dcj)) 

Si = —tti sin 9 + hi cos 9 cos </> + Cj cos 9 sin (j) , ti = —hi sin + Ci cos (j) . 



The two conditions to obtain the central configurations are there- 
fore 



Which are two equations to determine the two angles 9 and 

Starting from an initial value for the two angles in the interior of 
the spherical triangle allowed by the order of masses in the line, we 
iterate looking for the root of one angle assuming the given value for 
the other alternatively in the two previous conditions and then one 
iterates until desired convergence. 

Choosing the values of the masses as mi = 20, m2 = 13, ms = 
7, m4 = 6,1 have been able to compute the twelve different central 
configurations which are grouped in the table and represented his po- 
sition in the hemisphere of angles 9,(j) hy its stereographic projection 
in the figure 2. 

The four particles moving in a central configuration is illustrated 
by the last figure. The four masses move in similar conies that in 
figure 3 were selected ellipses of eccentricity 0.7. The orbits of the 
ellipses comes from integration of equations (28) and (29) as 



where e is the eccentricity and a is the major semi-axis of the ellipse. 
The constant of integration ipo measure the orientation of the orbits, 
and in figure 3 was chosen equal to 27r/5. The explicit coordinates for 



(33) 




(34) 



R = 




(35) 



1 — e cos(^ — ipo) 
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the four orbits were computed as function of tp from 
/ xi X2 X3 2^4 \ _ a(l — e^) / cosV' 

\ yi y2 ys yA j 1 - ecos(V' - V'o) V ^^"^ 

(36) 

The orbits are scaled from the elhpse R, ip by the four coordinates 
of the central configuration tabulated for the particular permutation 
13, 7, 20, 6. 
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